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We consider the persistence probability, the occupation-time distribution and the distribution of 
the number of zero crossings for discrete or (equivalently) discretely sampled Gaussian Stationary 
Processes (GSPs) of zero mean. We first consider the Ornstein-Uhlenbeck process, finding expres- 
sions for the mean and variance of the number of crossings and the 'partial survival' probability. We 
then elaborate on the correlator expansion developed in an earlier paper (G. C. M. A. Ehrhardt and 
A. J. Bray, Phys. Rev. Lett. 88, 070602 (2001)) to calculate discretely sampled persistence exponents 
C"j of GSPs of known correlator by means of a series expansion in the correlator. We apply this method 

O ' to the processes d n x/dt n — rj(t) with n > 3, incorporating an extrapolation of the series to the limit 

of continuous sampling. We then extend the correlator method to calculate the occupation-time 
and crossing-number distributions, as well as their partial-survival distributions and the means and 
variances of the occupation time and number of crossings. We apply these general methods to the 
d n x/dt n = rj(t) processes for n = 1 (random walk), n = 2 (random acceleration) and larger n, and to 
simple diffusion from random initial conditions in 1-3 dimensions. The results for discrete sampling 
are extrapolated to the continuum limit where possible. 



I. INTRODUCTION 



-a 

Stochastic processes driven by Gaussian white noise have a wide range of applications in the physical sciences and 
• beyond, ranging from Brownian motion to options pricing. Here we focus on two basic properties of a stochastic 
Gaussian time series: the number of crossings of the mean value of the series, and the fraction of time for which the 
series is above its mean value. The former, termed the crossing number, has long been of interest to engineers and 
mathematicians [1—3], and more recently to physicists [4,5]. The latter, termed the occupation time, has also been 
, studied by mathematicians for a long time [6-8] for both Gaussian and non-Gaussian stochastic processes and has 
recently seen a revival in the physics community in the context of nonequilibrium systems [9-13]. The occupation- 
t-H , time distribution for a stochastic process is also important due to its potential applications in a variety of physical 
systems which include optical imaging [14], analysis of the morphology of growing surfaces [15], analysis of temperature 
fluctuations in weather records [16], in disordered systems [17] and also due to the connection between the occupation 
time in certain discrete sequences and spin- glass models [18]. 

Of particular interest is a limiting case of these two properties, namely the probability, termed the persistence 
probability [19], that the time series is always above its mean value up to time T. The latter, for stationary Gaussian 

S time series, typically decays as exp(— 0T) for T large, with the exponent in general taking a nontrivial value. 
For continuous processes X(T), the Independent Interval Approximation (IIA) [20,21] may be used to calculate 
approximately (for continuous sampling) the asymptotic (i.e. large-T) forms of some of the probability distributions 
above. This approach, which makes the (generally invalid) assumption that the time intervals between zero-crossings 
are statistically independent, is surprisingly accurate in many cases. However, the IIA involves an uncontrolled 
approximation, which cannot be improved upon in general, and whose numerical accuracy is hard to estimate. Until 
now, the IIA has been the only general analytical technique available. In the absence of exact general results, 
calculations of probability distributions exist only for certain specific processes, although a short-time expansion for 
a general process has been developed [22] . 

Both continuous and intrinsically discrete time-series can be studied. In this paper we consider discrete-time 
sampling of an underlying continuous Gaussian stationary process (GSP), X(T), with zero mean, unit variance, and 
known correlator C{T) = (X(T)X(O)). We sample this process every time step AT, and study the discrete-time series 
X(i/XT). This is, of course, completely equivalent to studying a discrete process with the same correlators C(jAT). In 
[23,24] we have studied the persistence of a discretely sampled random walk and of a randomly accelerated particle, 
both of which can be mapped to a GSP by a change of variables. In [25] we calculated, using a series expansion in the 
correlator C(jAT), the persistence probability for an arbitrary discretely-sampled GSP. By extrapolating to the limit 
in which the time between samplings tends to zero, we obtained results for several continuum processes. Thus the 
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results developed for discretely sampled processes may, for sufficiently smooth processes, be extended to give results 
for continuous-time processes. 

Before going further, we illustrate the main ideas by considering the simple example of a stochastic process, the 
continuous-time random walk described by the Langevin equation 

x = £(t) (1) 

where is Gaussian white noise, (£(t)) = 0, and (£{t)£(t')) — 2DS(t — t'). The probability that x > up to time t 
decays as t~ e where the 'persistence exponent' is 9 = 1/2. This process is not stationary since its correlator, 

C{t u t 2 ) = ((x(h) - (x))(x(t 2 ) - (x))) = 2Dmin(ti,t 2 ) (2) 

does not depend only on the time difference, |ti — *2 j - Note that, since the process is Gaussian, it is completely 
specified by its correlator and mean. We can map the random walk onto a stationary process by a change of variables; 
we change to logarithmic time, T = hit, and to a normalized process X(T) via 

x(t) - (x(t)) 
y]{x{ff) (x(t)f 

obtaining the equation 

where r)(T) is again a Gaussian white noise, and X(T) has zero mean. The process X(T) is the Ornstein 
process. It is stationary, with correlator 

C{T U T 2 ) = (X(Ti)X(T 2 )) = exp(-|T 1 - T 2 |/2) . (5) 

An equivalent mapping can be made for any non-stationary Gaussian processes for which the correlator has the form 
C{ti,t2) — tfg{ti/t 2 )- Thus although here we only attempt to analyze stationary processes, the results are more 
widely applicable. Note that the exponential decay, exp(— 9T), is equivalent to the power-law decay, t~ e , hence the 
terminology 'persistence exponent' for 9. 

The occupation time is the number s of positive (say) values obtained from the n measurements of X(T). Let R n ^ s 
be the probability distribution of s for given n and r = s/n be the fraction of measurements that are positive. In 
the limit n — ► oo, s — ► oo, with r = s/n fixed, R n , s has the asymptotic form R n . s ~ [p( r )]" = cx p[ — Qd{t)T], where 
T = nXT and 0o(r) = — ln[p(r)]/AT. Here #d(1) = 9d(0) (all or none of the measurements positive) is the usual 
discrete persistence exponent introduced in [23]. In a similar way we can define P n , m to be the probability of observing 
m zero-crossings in n measurements. If now r = m/n and we take the limit n — > oo, m — > oo, holding r = fixed, 
we find P n , m ~ [p( r )]™ = ex P[ — Qd{t)T\. Here 6d(0) ( n0 crossings) corresponds to the usual discrete persistence. 
Although we use the same symbols p{r) and 9n(r) for the occupation-time and crossing problems, it should be clear 
from the context which problem we are referring to. 

In this paper we extend the technique of [25] to calculate the exponents 6d(t), or equivalently the functions p(r), 
for the occupation-time distribution and the distribution of crossings for arbitrary discrete or discretely-sampled 
Gaussian stationary processes. The technique gives the exponents as a series expansion in the correlators C(j/XT) 
up to C(10/XT) and C(XT) 10 , i.e. 10th order. For the calculation of the persistence exponent we work to 14th order. 
The results work well for C(j'AT) small, i.e. the time between samplings large compared to the correlation time of 
the stationary process. For certain processes we are able to extrapolate the series to the limit XT — > 0, thus obtaining 
values of the continuum exponents that compare favourably with those predicted by the IIA when measured against 
exact or numerical results. 

The layout of this paper is as follows: In Part I we consider the Ornstein-Uhlenbeck process introduced above, 
this being perhaps the simplest of GSPs. By extending the 'matrix method' developed in [23] we find the mean 
and variance of the distribution of zero-crossings as a perturbation expansion to high accuracy. We also use another 
method to find the same results. The case of an unstable potential (obtained by changing the sign of the drift term) is 
also considered. We extend the concept of partial survival [4] to discrete sampling, calculating F n (p), the probability 
of surviving to n samplings if each detected zero-crossing is survived with probability p. 

In Part II we apply and extend the correlator expansion method of calculating the persistence exponents, occupation- 
time distribution and crossing distribution of an arbitrary discretely-sampled GSP. The general method is then applied 
to some specific examples of interest, and the results extrapolated to the limit of continuum sampling where possible. 
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In section VI we introduce the correlator expansion first developed to 14th order in [25] as a method of calculating 
discretely-sampled persistence exponents. We explain this technique more fully including the extrapolation to the 
continuum limit using constrained Pade approximants, which allows rather accurate calculation of the standard 
persistence exponents. In [25] this technique was applied to the case of the random acceleration process and also to 
diffusion from random initial conditions in 1-3 dimensions. Here we also apply it to the class of processes d n x/dt n = 
rj(t) where 77(f) is Gaussian white noise. The n = 1, 2 cases are the random walk and random acceleration problems 
already studied. Here we calculate the persistence exponents for larger values of n and show numerically that 
Q n — #oo oc 1/n for n large. The results are compared to the predictions of the Independent Interval Approximation 
(IIA). 

In Sections VII and VIII we extend the correlator expansion to calculate the occupation-time distribution, i? n?s , 
this being the probability of s positive measurements in n samplings, to 10th order in the correlator. This gives, in 
particular, the variance of the occupation-time distribution which we calculate in a more straightforward way as a 
check. We define a partial-survival occupation probability, P n (p), as the probability of surviving to n samplings if each 
positive sampling is survived with probability p. This is also the generating function for i?„ jS . We find P n (p) to 14th 
order in the correlator. We apply the results to the following five GSPs: the random walk, where the results of Part I 
are used as a check of the method, the random acceleration problem and diffusion from random initial conditions in 
1, 2 and 3 dimensions. Extrapolations to the continuum are included. 

In Sections IX, X and XI we further extend the correlator expansion to calculate P n ,m, the probability of m detected 
crossings in n samplings. This is found to 10th order in the correlator and also enables us to calculate the partial- 
survival probability and the moments of the crossing distribution. We apply the result to the five GSPs of section 
VII, and also to the d n x/dt n = 77(f) processes for n > 2 and to an intrinsically discrete process for which the exact 
results are known [5]. Extrapolations to the continuum are included and we compare continuum results with the 
IIA and also, for the random acceleration partial-survival problem and the intrinsically discrete process, to the exact 
solutions. In Section X we show the result for the mean number of detected crossings. We also derive the variance as 
a series expansion in the correlator, the coefficients of which agree with those of the previous section. We conclude 
with brief summary of the results. 



Part I 

The Ornstein-Uhlenbeck Process 

In Part I we will study the detected crossings of the Ornstein-Uhlenbeck process. Beside being of interest in its own 
right, this will illustrate some of the methods used later and also provide several checks on the correlator expansion of 
Part II. Furthermore, here we are able to calculate probability distributions starting at a certain position X , rather 
than just the long-time or stationary state distributions. 

Consider the stationary Gaussian Markov process evolving via the Langevin equation, 

^L = -^X + V (T), (6) 

where r/(T) is a white noise with mean zero and correlator, (r)(T)i](T')) = 2DS(T—T'). This is the Ornstein-Uhlenbeck 
process whose persistence exponent for discrete sampling was calculated in [23]. 
Integrating equation (6), we get 

X(T) = X e-» T + e-» T f r){T x )e^dT u (7) 
Jo 

where Xq = X(T = 0). Let T = nAT. Then the mean position (X n ) after n steps starting initially at Xq is given 
from equation (7) by, 

(X n ) = X e-» nAr = X a n , (8) 

where a — e _AlAr . Similarly one finds that the correlation function is given by, 

([X n - (X n )][X m - (X m )}) = D>(a} n - m \ - a"+ m ), (9) 

where D' = D j '/i. Thus in the stationary state, n — ► oo, m — > oo with n — m fixed, this process has mean zero and a 
correlator C{T U T 2 ) = (V(Ti)V(T 2 )) given by 

C(T U T 2 ) = <7(|T 2 -Ti|) = D'e-^~ T ^. (10) 
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II. DISCRETE BACKWARD FOKKER-PLANCK EQUATION 



Let Q n (m 7 X) be the probability that starting at X at T = 0, the process, when sampled only at the discrete 
points, changes sign m times within n discrete steps. Note that the probability P n ^ m of observing m sign changes 
in n measurements (as defined in the introduction) can be simply obtained from Q n (m,X) via the relation, P n;m = 
JTinfty Qn(iri 7 X)p (X)dX where po(X) is the initial distribution of X which we will take as the stationary distribution 
of X. One can write down a recursion relation for Q n (m, X), valid for X > 0, by noting that at the first step either 
the process changes sign or it does not. 

/>oo />0 

Q n+1 {m,X)= Q n {m,Y)G{Y,XT\X,0)dY + Q n {m-l,Y)G(Y,XT\X,0)dY, (11) 

JO J-oo 



where G(Y, /XT\X, 0) is the probability of going from X to Y in a time AT, given by 

I (Y-aX) 2 



G(Y,£T\X, 0) 



v /2ttL>'(1 - a 2 ) 



g 2D'(1- 



(12) 



Using the rescaled variables x = Xj y/D'(l — a 2 ) and y = Y j D'{\ — a 2 ), and making use of the symmetry 
Q n (m, -x) = Q n {m, x) we get, 



Q n+1 (m,x) 



Q n (m,y)e' 



-(y-ax) 2 /2 



+ Q n (m-l,y)e-(y +a ^ 2 / 2 



dy 



(13) 



Note that for m = this reduces to the persistence problem studied in [23], whilst for m — n we have the 'alternating 
persistence' problem [23]. 

We define the generating function, F n (p,x) — ~Y^ n= r > Qn{ l m,x)p rn ' . The generating function has a physical inter- 
pretation: If one considers that with every detected change of sign a particle survives with probability p (partial 
survival), then F n (p, x) is precisely the survival probability of the particle. Note that for p = 1 this probability is 1 
whilst for p — we recover the usual discrete persistence. 

Multiplying equation (13) by p m and summing over m we get, 



F n+1 (p, x) 



1 f°° 

75? i, F " lp - V) 



e -(y-ax) 2 /2 +pe -(y+ax) 2 /2 



dy. 



(14) 



Let E n {x) = X)m=o m Qn(m, x) = dF n {p,x)/dp\ p -\ denote the expected number of sign changes up to n steps 
starting at x at T = 0. Taking the derivative with respect to p and putting p = 1 in equation (14), it follows that 
E n (x) satisfies the recursion relation, 



E n+1 (x) = 



E n (y) 



e -(y-ax) 2 /2 + e -(y+ax) 2 /2 



dy + -erfc 



ax 
V2j' 



(15) 



where erfc(x) is the standard complimentary error function and the recursion in equation (15) starts with the initial 
condition, E (x) = 0. Also note that G n (x) = J2m=o m { m ~ ^)Qn(m,x) = d 2 F n (p, x)/dp 2 \ p=1 satisfies the recursion, 



G n +i(x) 



/2tt 



G n (y) 



e -(y-ax) 2 /2 + e -(y+axf/2 



dy + 



dyE n (y)e~ 



-(y+ax) 2 /2 



(16) 



with the initial condition Gq{x) = 0, and where E n [x) is given by the solution of equation (15). In order to calculate 
the average number of crossings and the variance around this average, we need to solve the two integral equations 
(15) and (16). 

If we can obtain the solution of E n (x) from equation (15), then we need to average over the distribution of the 
initial position x to obtain (m) n — E n (x)po(X)dX where Po(X) is the initial distribution of the position X 



and x = X/y/D'(l 



For a < 1, i.e. p, > 0, if we choose po(X) to be the stationary distribution of the 



process, po(X) = = e x l 2D , which corresponds to waiting an infinitely large number of steps before starting 

the measurements, then (m) n can be computed exactly from equation (15). Multiplying equation (15) by this Po{X) 
and integrating over X we get 



( m )n+l = ( m )n + A > 



(17) 
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where A = | — -^sin 1 (a). Using (m) = 0, we get exactly, 

(m>„ = n\ (18) 

which is the random walk case of the general result (m) n /n = 1/2 - (1/tt) sin -1 C(AT) [26]. Note that for AT -> 
this reduces to Rice's formula, while for AT — > oo it reduces to n/2 as expected since the values of X at the discrete 
points become statistically independent. Note that for a > 1, there is no stationary distribution as this corresponds 
to an unstable potential. 

Note however that to compute g n — (m(m — l)) n by a similar method, we need to know the full function E n (x), 
that is we need to solve the full integral equation (15). Indeed for a < 1, if we choose to average over the stationary 
distribution, po(X) = -p===e~ x l 2D , then by multiplying both sides of equation (16) by po(X) and integrating over 
x from — oo to oo, we get after straightforward algebra 

9n + i = 9n + \[^^[ En(y) e^^^erfc dy. (19) 

Hence the variance of the number of crossings, a 2 = (m 2 ) n — (m) 2 n is given by, 

crl = g n + n\-n 2 \ 2 , (20) 

where g n is given by the solution of the recursion equation (19) and A = ^ — isin _1 (a). To determine g n from equation 
(19), we need to know the full function E n (x). 

In the next section, we show that there is an alternative way to derive an expression for E n (x) without solving the 
integral equation (15). 



III. ALTERNATIVE DERIVATION OF E N (X) 



An alternative derivation of E n (X) can be obtained by noting the evident relation, 

E n+1 (X) - E n (X) = \[l~ (sign(X n )sign(X n+1 ))] 

= (6(X n )) + (6(X n+1 )) - 2(6(X n )6(X n+1 )), 



(21) 



where 9{x) is the standard theta function. Noting that X n and X n+ \ are Gaussian variables, we can compute the 
right hand side of equation (21) exactly. For this we first need the joint distribution P[X n , X n+1 ] of X n and X n+ i 
which involves the correlation matrix C n>n +i given by 



+ 1 



D' 



l-a 2n a(l-a 2n ) 
a(l-a 2n ) l-a 2n + 2 



where we have used equation (9) for the matrix elements. Using this joint distribution and carrying out the Gaussian 
integrations, we get after lengthy but straightforward algebra the final expression of the right hand side of equation 
(21) as, 



E n+1 (X) - E n (X) = i e rfc(-A n X) + ^erfc (-A n+1 X) 



_l r 

J-A n+1 X 



dye y erfc 



1 -a 



2n+2 



I -a 2 



-A n X 



l-a 2n 

- —y 



l -a 2 



(22) 



where A n = a n /y / 2D'(l — a 2n ). Changing to rescaled variable, x = X/yjD'{l — a 2 ) and using the initial condition 
E (x) = 0, we get from equation (22), 



_^ n— 1 



m=0 . 



erfc 



a m x 1 — a 2 



V2 V 1 - a 2 



+ erfc 



V2 



I -a 2 



1 - a 2m + 2 



2 r°° 

\pH J_ a m + 1 x / 1-a 2 



■ dye v erfc 



1 - a 2m + 2 1-a 2 a 



1-a 2 V 1 - a 2m 



1 - a 2m 
7^ 



-y 



(23) 
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This can be solved numerically, figure (1) shows E n (x) — nX for a = 1/2 and n = 1000. 
0.5 
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FIG. 1. Plot of u(x) = E n (x) - nX for a = 1/2 and 
n — 1000, calculated using equation (23). u(x) is symmet- 
ric about x = . Also shown is u(x) for a = 1/2 derived 
using the matrix method (equation(30)), the two curves are 
indistinguishable. 

In the next section we discuss the exact asymptotic properties of the mean and the variance of the number of sign 
changes for the two cases: (A) Stable potential, \i > 0, i.e < a < 1 and (B) unstable potential, \i < 0, i.e. a > 1. 

IV. ASYMPTOTICS OF THE MEAN AND VARIANCE OF THE NUMBER OF SIGN CHANGES 

A. Stable potential: fj, > 



In this case a = e~ MAT < 1 and also D' > 0. Thus as n — > oo, A n = a n / y / 2D'(l — a 2 ) — > 0. Then from equation 
(22), we find E n+ \{X) — E n (X) — > A with A = ^ — isin~ 1 (a) indicating that in the n — > oo limit, E n (X) becomes 
independent of X (as long as as X < a n ) and to leading order for large n, 



E n (x) = nX, 



(24) 



in agreement with equation (18). Indeed the expression for E n (x) in equation (23) is a solution of the integral equation 
(15). For finite n, one can write E n (x) = nX + u n (x). An explicit expression for u n {x) can be obtained from that of 
E n (x) in equation (23). Using the explicit value of A and after a few steps of algebra we get, 



n-l 



ln{x) = - S m (x 



m=0 



s m (x) = erfc ( -^JJ-^j + crfc ^-i-lf^l-^ - 2 + -| ^ d ye ^ 2 erfc ( 



V2 



-ay 



2 



» r° 



V2 



+ 1 / l-a 2 « 



■dye y erfc 



l_ a 2m + 2 



l-a : 



2m+2 



l-o 2 a" 



1 - a 2 V 1 - a 2m 



a; — a 



l-a 2r 
1-a 2 



(25) 



Besides, using E n (x) = nX + u n (x) in the integral equation (15), we find that u n (x) also satisfies the following integral 
equation, 



u n+ i(x) 



1 f°° 
= ~7k= / Un 

V2^y 



(y) 



e -(y-ax) 2 /2 + e -(y+aa;) 2 /2 



(26) 



G 



with A = 5 



isin-Va). 

7T V / 

Now the leading term in E n (x) is nA and is independent of x as long as x < a~ n . This upper cut-off tends to oo as 
n — > oo since a < 1. The a; dependence of E n (x) appears only in the subleading term u n (x). Now, as n — ► oo, u n (x) 
tends to a stationary solution independent of n (as long asi< a - ") and is given by the fixed point solution u{x) of 
the integral equation (26), 



u(x) 



1 



e -(y-ax) 2 /2 + e -(t/+a a; ) 2 /2 



'/</ + ^erfc ( 



A. 



We can find u n (a;) perturbatively by expanding the axy term in the exponentials in equation (27) to get, 
1 



u(x) — 



(ax) m e-^ I 2 f°° . , 2, 1 

^ : / u(y)e y /2 2y m dy + -erfc 

Z7T , . mi J i 

m—0 (meven) 



E 



ax 
V2 



A. 



(27) 



(28) 



Defining 



gives 



-/ 

Jo 



u(x) 



a l ' 2 x l 



- x2 ' 2 dx 



u(x) 



E 

m—0 (meven) 



jirt^mg- (ax) 2 



"ml 



m+ 2 



A. 



(29) 



(30) 



Multiplying both sides by a l / 2 x l e x2 / 2 /W. and integrating over positive x gives, 

Ji 



ii= M ^ J " 

m—0 (meven) 



(31) 



where 



M lm = 



1 



2a 



2ira V 1 + ° 2 



(i+ m +l)/2 



r[(i + m+l)/2] 



m even, otherwise 



-erfc 
2 



Thus 



I = (1 M) _1 J 



(32) 
(33) 

(34) 



This perturbative expansion in a is the matrix method. Although M\ m is an infinite array, the elements decrease 
rapidly with increasing I and m since each increment of I and m gives a higher power of a (a < 1). We can solve 
this numerically. Figure 1 shows u(x) for a = 1/2. Alternatively, using mathematica, we can solve equation (31) 
iteratively. At each iteration we obtain a new I whose elements are a series in a up to our required order. Convergence 
is rapid. This can only be done with a smaller matrix than the numerical method, but gives a result for general a. 
Note that we could equivalently have done this perturbative expansion on equation (15) to get E(x) and subtracted 
the nA. 

Alternatively, we can directly take the n — > oo limit of the expression of u n {x) in equation (25) to get, 



1 x 

*0) = 9 E S m{x), 



(35) 



m=0 



where s m (x) is given by equation (25). This has been done numerically (for n large but finite) and found to agree 
with the matrix method. 

Let us first compute the asymptotic properties of the fixed point solution u(x). Consider first the limit x — > 0. 
Putting x = in equation (25) and carrying out the integrations we get, after some algebra, 



7 



1 CO 



sm 



m=0 



7^ 



,2m+2 



(l - v 7 ! - a 2m ) 



(36) 



For example, for a = 1/2, we get u(0) — 0.19160374 . . . which agrees very well with the result obtained from the direct 
numerical integration of equation (26) in the large n limit. Consider now the other limit x — > oo. By making the 
change of variable, y — ax = \/2z in the integral equation (27), we get to leading order for large x (where the lower 
limit of the first integration — > — oo), 



u(x) w u(ax) — A. 



The solution of this equation is given by, 



u(x) 



In a 



In a;. 



(37) 



(38) 



Note that A/In a < for a < 1 and hence u(x) goes to — oo logarithmically as x — > oo. This is consistent with the fact 
that E n (x) w n\ + Aj^j <~ jj^ln(ia n ) — > as x — > a - " as it should evidently from the direct expression of E n (x) in 
equation (23). 

In figure 1, we plot u(x) obtained from numerically evaluating the sum in equation (25) and also the result obtained 
using the matrix method. The results agree to within numerical precision. Also, they agree with the asymptotic 
results in the large and small x limits. Note that the x — > oo limit was not attainable by the matrix or summation 
methods because in both cases the evaluation is done to a finite order or finite n. 

Once we know u n {x), then we can determine the variance cx 2 from equation (20). Substituting E n (x) = n\ + u n (x) 
in equation (19) and using the fact that u n (x) tends to the fixed point solution u n (x) — > u(x) for large n, we get 



9n+i = 9n + 2nA 2 + (3, 



where 



W Z£l J™ U (y) e -^i--)y»/»erfc {^j dy. 



(39) 



(40) 



One can easily solve the recursion equation (39) exactly using g = and we get, 

g n - \ 2 n 2 + (/3- A>, 

where [3 is given by equation (40). Substituting this expression for g n in equation (20), we finally get the required 
exact expression of the variance for large n, 



(41) 



(A - A 2 + 0) n. 



(42) 



Thus the variance can be exactly determined once we know the function u(x) and thereby [3 from equation (40). 
Using the exact expression of u(x) from equation (35), we have, in principle, an exact result for j3 and hence for <r 2 . 
Substituting the u(x) derived from the matrix method into equation (40) gives (3 and hence cr 2 . cr 2 /n is plotted as a 
function of a for D' = 1 in figure (2). 
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FIG. 2. Plot of ol/n against a = e~ MAT with D' = 1. Note 
that (Tn/n — » oo for a — » 1, and the series has not converged 
for large a. 



B. Unstable potential: /i < 

In the previous subsection, we have seen that for a < 1, the function E n (x) behaves asymptotically for large n 
as E n {x) = nX + u(x) where u(x) is given cither by the exact expression in equation (35) or equivalently by the 
solution of the integral equation (27). For the unstable potential /x < 0, i.e. a > 1, the number of crossings will be 
finite and so E n (x) approaches a steady state as n — > oo. This is most easily seen from equation (22). For a > 1, 

oo (note that D' = D//i < 0). Taking this limit in equation (22), we see that 

oo. This steady state E(x) is given by 



A n = a n /^2D'{\ - a 2 ) -> oo as n 

E n+ i(x) — E n (x) — > for all x as n — > oo, indicating E n (x) — > as n 

the fixed point solution of the integral equation (15) with a > 1, 



1 



e -(f— ) 2 /2 + e - 



(j/+aa;) 2 /2 



1 r ( ax 
-errc 

2 VV2 



i?(x) can be found using the matrix method in the same way as before but with J; replaced by J/ where 

■ /,_ W 

_E(x) is shown in figure (3) for the case a = 2. 



x l e 



-^erfcf^ 



V2 



(43) 



(44) 




FIG. 3. The expected number of detected crossings E(x) 
for the case a = 2 (unstable potential). Shown are the re- 
sults of the matrix method and also numerical evaluation of 
equation (45) which are indistinguishable. 



Alternatively, E(x) is also given by taking the n — > oo limit of the exact expression in equation (23) (with a > 1), 

a m x 



1 oo r / 
(x) = - erfc 

m=0 L V 



erfc 



a m+1 a; 



V2 



1 



1 - a 2m + 2 



erfc 



'1 - a 2m + 2 / 1 



1 - a 2 Vl-o 2ra V2 



.x — a 



7 2m 



1-a 2 



(45) 



This is shown in figure (3). 

We have calculated the variance of the number of detected crossings and also the expected number of detected 
crossings starting at position x. These calculations have been carried out by two independent methods and the results 
agree with other. 
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V. PARTIAL SURVIVAL 



The partial-survival probability, F n (p, x) , is the probability of surviving beyond the nth sampling having started 
at x if each detected crossing of the origin is survived with probability p. Thus, 

oo 

F n (p,x)= J2Qn(m,x)p m (46) 

and as stated before F n (p, x) is the generating function for Q n {m, x). F n (p 7 x) satisfies the integral equation (14). We 
expect that for large n, F n (p,x) = [p p (a)] n F(x) where p p (a) — e~ e ^^ and 8{p) is the discrete persistence exponent. 
Substituting this into equation (14), we get an eigenvalue equation for F(x), 

p p (a)F(x) = ^ F{y) e- {y - ax ^l 2 +pe-^ +ax '> 2 ' 2 dy. (47) 

V27T Jo L J 

The largest eigenvalue p p {a) and the corresponding eigenfunction can then be determined either by the matrix method 
or by the variational method as in our previous paper [23]. Using the matrix method, we get, 

Pp (a)F(x) = e -^=- £ i-a^l + (-l)»/ m (48) 

v 771 — 

OO 

p p (a)I l =Y t {G lm +pG , lm )I m (49) 

m=0 

where 

m roc 

l m = — x \ dyy m e-y / 2 F(y) (50) 
ml J Q 

g „ ■ r(*\^ w m±z+m (51) 



8tt \l + a 2 J l\m\ 

and G' lm (a) = Gi m (—a). In fact, G/ m is the matrix used for calculating the discrete persistence exponent [23], whilst 
G' lm gives alternating persistence. This is to be expected, since for p = we just have ordinary persistence and for 
p> 1 we would expect the paths which cross between every sampling (alternating persistence) to dominate. In the 
same way, one may calculate p p (a) for the discretely sampled random acceleration problem studied in [24], whose 
stationary process is given by 

X + (a + p)X + apX = 7](T) (52) 

where i){T) is Gaussian white noise with mean zero and correlator < rj{T)rj{T') >= 2a[3{alpha + [3)5(T — T"), and 
a = 1/2, /3 = 3/2 for the random acceleration problem, although other values of a, [3 can be considered. We get 

oo 

p p (a)Iij = ( H ijki + pH'ijki) hi (53) 

k,l=0 

where Hijki and Jjj are given in [24]. Again, Hijki is the matrix used to find the discrete persistence exponent and 
H'ijki( a ) = Hijki(-a) gives the alternating persistence exponent. We find p p (a) numerically and also as a power series 
in a for the two processes given above. The results are shown in figures (4,5). Also the eigenfunction F(x) for the 
Ornstcin-Uhlcnbcck Process may be found by substituting the eigenvector corresponding to the largest eigenvalue 
into equation (48). Results for a = 0.5 with p = 0, 0.5, 1 are shown in figure (6). 
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0.5 ^ 1 1 1 1 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

exp(-nAT) 

FIG. 4. Plot of the random walk discrete persistence eigen- 
value p P (a) for partial survival against a = e~ pAT for values 
of the survival probability p from (normal discrete persis- 
tence, lowest curve) to 1 (guaranteed to survive so pi (a) = 1, 
top curve) in steps of 0.1. The curves are the raw series in a 
to order a 50 . Note that for all the curves, p p — > 1 for a — > 1 
since a walker will always survive for a time AT when AT — > 0. 
Since the series in a are finite, they do not quite converge to 
1 in this limit 



o. 




0.2 0.4 0.6 0.8 1 



exp(-AT/2) 

FIG. 5. Plot of the random acceleration discrete persis- 
tence eigenvalue p p (a) for partial survival against a = e _AT ^ 2 
for values of p from (normal discrete persistence, lowest 
curve) to 1 (guaranteed to survive so pi(a) — 1, top curve) in 
steps of 0.1. The curves are Pades of the raw series in a to 
order o . 



5 10 15 20 

FIG. 6. The eigenfunctions F(x) of equation (47) for 

a — 0.5 with p = 1 (lowest curve), p = 0.5 (middle curve) 
and p = top curve. The eigenfunctions are defined only up 
to an arbitrary prefactor which has been chosen here so that 
F(0) = 1. Since the eigenfunctions are series in x, for large 
x they do not converge to the correct solution. This can be 
clearly seen for the p = 1 case where, since the walker is guar- 
anteed to survive, F(x) is a constant everywhere whereas the 
plot is not constant for x > 10. For the p = case, F(x) ~ x v 
with v = hip/ In a. 



So far we have found the mean and variance and also the partial-survival probability for the Ornstein-Uhlenbeck 
Process, a simple Gaussian Markov problem. We did this by using the propagator, P(Y, /XT\X 7 0) (equation(12)). 
We also showed how the partial-survival probability of other GSPs of known propagator can be found by using the 
perturbative matrix method, and we illustrated this for the random acceleration problem. However, the methods used 
become progressively harder as the number of variables in the problem increases. The Ornstein-Uhlenbeck Process 
had only the position X, the random acceleration problem had X and V, while the persistence problem for the 
diffusion equation cannot be expressed in terms of a propagator with a finite number of variables. In the remainder of 
this paper we will use a different approach based on the correlator of the process C (T) . This removes the difficulties 
mentioned above and can, furthermore, be applied to any GSP of known correlator. The results for low-dimensional 
problems obtained above are slightly more accurate than those given by the correlator method because they can be 
calculated to higher order. They can be used as powerful checks of the accuracy of the correlator method. 
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Part II 

The Correlator Expansion 

The correlator expansion was initially used to calculate the discrete persistence exponent of an arbitrary Gaussian 
Stationary Process (GSP) and also, through extrapolation to the continuum, the continuum persistence exponent 
[25]. Here we will extend the method to calculate the occupation-time distribution and the distribution of crossings. 
As these calculations will require some explanation, we take this opportunity to describe the correlator expansion in 
full. 



VI. AN INTRODUCTION TO THE CORRELATOR EXPANSION 



The expansion starts from the following identity for P n , the probability of no detected crossings in n samplings: 

P n = (j[0[X(iZT)}j (54) 

where @(X) is the Heaviside step function and the expectation value is taken in the stationary state. One may write 
Q[X(iAT)] = (1 + <7j)/2, where Oi = signLY(iAT)], and expand the product to give, 

\ 

+ X! (piOjOkOi) + • • • J (55) 

l=i<j<k<l J 

where the terms with odd numbers of as vanish since the process is symmetric under X — ► —X (and therefore under 
a — > —a. To evaluate the terms we use the representation 

1 [°° dz t Zl e tz > x > 

ai = — hm / — — — (56) 

nr e^o J_ QO (zi - ie){zi + le) 

Carrying out the required averages of the Gaussian process gives the correlation functions appearing in (55): 

(^i • • • <0 = / fl ( ^7 ) exp ( ~\ Za Ca P z f } ) ' ( 57 ) 




where C Q( 3 = (X[aAT] X[(3/XT\) = C(\a — /3|AT), and there is an implied summation over a and (3 from 1 to m. 
Notice that we have already taken the limit e — > in (57), with the understanding that all integrals are now principal 
part integrals. 

For the m = 2 case this integral can be done exactly by differentiating with respect to C\i and doing the two 
simple Gaussian integrals before integrating again with respect to C\i and imposing the boundary condition that 
(cj 1 o'i 2 ) = for C\2 — 0. This gives the well-known result (cr^a^) = (2/n) sin -1 C\i. For m > 4 this method becomes 
non-trivial. Instead, we choose to expand the exponential in equation (57) in powers of C a p (a ^ (3) leaving the terms 
with a = (3 unexpanded (noting that C aa = 1). This allows us to evaluate each correlation function of the as up to a 
given order in the correlators C a p. By symmetry, only terms which generate odd powers of every z a in the expansion 
of the exponential (to give even powers overall in the integrand, through the factors 1/zi) give a non-zero integral. 
This suggests a simple diagrammatic representation for the terms in (55), as given by (57). On a one-dimensional 
lattice containing n sites, with lattice spacing AT, draw m vertices at the locations l\, I2, ■ ■ ■ , l m - Connect the vertices 
by lines in all possible ways (summing over these different possibilities) subject to the constraint that each vertex is 
connected to an odd number of lines. Associate a factor \/~2n(p— 2)!! (coming from evaluating the Gaussian integrals) 
with each vertex of order p, a factor (— Ci^Y /r\ with the r lines connecting site k to site lj, and an overall factor 
(7ri)~ m with the diagram. This suffices to evaluate the integrals in (57). Evaluating the summations in (55) involves 
enumerating all configurations of the vertices on the lattice for a given ordering of the points, and noting that the 
factor Ci i i j associated with a given line is equal to C(qAT), where q — \k — lj \ is the length of the line in units of AT. 

We choose to define C(AT) as 1 st order small and C(q/XT) as q th order small. Although somewhat arbitrary, this is 
a large AT expansion and most processes of physical interest have correlators C(qAT) which decrease exponentially 
for large argument, so our definition is consistent for large AT. For the random walk, the correlator is C(T) = e^^/ 2 
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and the definition is always valid. For other processes it is often possible to re-expand the correlator in terms of an 
exponential and work to a given order (in practice we can go up to 14th order) in this exponential. This can be 
done for, e.g., the random acceleration problem. The order of a diagram is then equal to the total length of its lines 
(measured in units of the lattice spacing AT). Thus to a given order k, we need only evaluate correlations functions 
with separations up to 2k AT. 

To illustrate this approach, we show in Figure 7 all the topologically distinct diagrams contributing to {<Ji<Jj(JkVi) 
up to 4 th order. The first diagram, when enumerated on the lattice, will be 2 nd order or greater, while the remaining 
five will be A th order or greater. 



FIG. 7. All topologically distinct contributions to 
(oiOj<Jk<Ji) up to fourth order. Note that the 1st, 5th and 
6th diagrams are disconnected, whilst the others (including 
diagrams 2 and 3) are connected (due to the constraint that 
the order of the points k,l must be unchanged). 



In Figure 8 are shown the enumerations of two of the basic diagrams of figure 7 together with their embedding 
factors (the number of ways they can be placed on the lattice), up to 5th order. 




+ (n-4) • <J)» • • + (n- 4 ) < 

Enumeration of two of the diagrams from figure 
7, with embedding factors, up to fifth order. Large dots are 
vertices, small dots intermediate sites. The final diagram here 
gives a contribution (n - 4)(2/7r 2 )C(2Ar)C(Ar) 3 to P n . 



Thus the calculation of P n proceeds in 3 stages: 

1: all the basic diagrams up to the required order are found, for example the 4- vertex diagrams shown in figure 7. 
2: the basic diagrams are enumerated on the lattice, including the 'stretched' diagrams, as shown in figure 8. 
3: the appropriate factors are assigned to each diagram. 

The total number of enumerated diagrams increases roughly by a factor of 2 for each extra order. At 14th order 
there are 12 434 diagrams. The process was automated using MATHEMATICA. For calculating P n , finding the basic 
diagrams was the most challenging task in terms of computer time and memory. To achieve 14th order we used 
the fact that all disconnected diagrams can be constructed by combining two or more connected diagrams. At 14th 
order, diagrams with 12 or more vertices are disconnected, thus only connected diagrams with up to 10 vertices need 
be found. Furthermore, diagrams containing 2 or more lines connecting the same points can be constructed from 
diagrams with only or 1 lines connecting points by adding pairs of lines. The procedure adopted was as follows: 

For the 2 to 10-vertex diagrams: all possible diagrams up to 14th order with only or 1 lines connecting points 
and all vertices odd are constructed. To these diagrams pairs of lines are added in all possible ways up to 14th order. 
The connected diagrams are then selected and stored. 

For 2 to 28-vertex diagrams: all diagrams are constructed by combining the connected diagrams found above, for 
example, the 6- vertex diagrams are formed from 3 2-vertex diagrams, one 2-vertex and one 4- vertex diagrams, and one 
6-vertex diagram (with the appropriate permutation factors arising from the various ways of ordering the connected 
diagrams). For large vertex numbers, this produces significant savings over naively trying all diagrams (since the 
vast majority of diagrams do not satisfy the odd- vertex criterion). In this way, up to order k we need only find all 
connected diagrams with up to q vertices, where q = (2k + 4)/3 is even and we round down, whilst the diagrams go 
up to 2k vertices. 

Having found P n , we find p(= e~ eAr ) using the fact that P n ~ p n for large n. Thus p is formally obtained as 
p = linijj—xx, P n+ i/P n . However, since we started in the stationary state, the relation p = P n+ \/P n in fact holds for all 
n larger than the length of the longest diagram. Expanding the expression for P n+ \/P n as a Taylor series up to 14th 
order ina = e -^/ 2 gives us a series expansion for p(a) about AT = oo (a = 0). For the random walk, for example, 
C(AT) = e^^/ 2 — a for p = 1/2, substituting this into our expression for p gives us a series up to 14th order in a 
whose coefficients agree with those of the matrix method (to within the numerical error of the matrix method). 

For the usual random acceleration problem, C(AT) = (3e _AT / 2 — e~ 3AT / 2 )/2 = (3a — a 3 )/2. For this case our 
identification of C(14AT) as being of the same order as C(AT) 1A does not hold for all AT. However, if we only keep 
terms up to a 14 our expansion be exact up to order 14 in a. Whenever possible, this is what we will always do. Note 
that in this way we are now working strictly to 14th order, even though C(jAT) ^ C(ATy . 
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In [23] we introduced the concept of alternating persistence, with being the probability that is positive for 
odd i and negative for even i (or vice versa). We find pa by noting that, whereas before we required X\, X2, ■ ■ ■ X n > 0, 
we now require X\, —X^,X^, —X4 . . . X n > 0. Thus the calculation is as before except that C(qXT) — ► —C(qXT) for 
q odd. Making this minor change to the normal persistence result gives us the alternating persistence exponent. This 
way of accounting for sign changes between X if Xj will be used below to calculate the distribution of crossings. 

We have applied the correlator expansion to the random walk, x = rj(t) and random acceleration, x = rj(t) using 
the transformation to logarithmic time to generate the corresponding stationary processes. Furthermore we have 
studied diffusion from random initial conditions in 1-3 dimensions, d<j>/dt = V 2 0, where </>(x, i) is the diffusion field 
and the initial condition </>(x, 0) is delta-correlated Gaussian noise. We consider the persistence of at a single site, 
for example <j>(0, t). For this process the correlator is: 

C{T) = sech d/2 (T/2) (58) 

where d is the space dimension. As for the random acceleration, we define a = exp(— AT/2) for d = 2 and a = 
exp(— /XT /A) for d = 1, 3 and then expand the correlator in powers of a. For d = 1, 2 the lowest power of a is a 1 and 
so we expand up to a 14 whilst for d = 3 the lowest power is a 3 and so we expand up to a 42 . We also considered the 
processes d n x/dt n = i](t) for n > 3. In logarithmic time the correlators are [20]: 

C n (T) = (2 - l/n)e- T / 2 2 F 1 (l, 1 - n; 1 + n; e- T ) (59) 

where 2^1 is the standard hypergeometric function. The n = 1, 2 cases are the random walk and random acceleration, 
whilst the limit n — > 00 case reproduces the correlator for the d = 2 diffusion process mentioned above [20] . 

For all these problems we define a discrete persistence exponent, 60(a) = — \np(a)/XT = mp(a)/2 In a, and plot 
9d{o) against a for < a < 1, i.e. 00 > XT > 0. Since we are plotting finite series in powers a, they do not converge 
for a — > 1. This problem is exacerbated by the 1/lna term in the definition of 9 u(a), which causes #d(1) to blow up 
unless p(l) = 1. To make p and hence 9d{cl) more accurate for a close to 1 we extrapolate #o(a) to the continuum. To 
do this we use the technique of Pade approximants borrowed from the field of series expansions for critical phenomena 
[27]. The Pade approximant involves replacing the 14th order series in a with a fraction whose numerator and 
denominator are series in a. The sum of the order of these two series is 14 and the coefficients are chosen so that 
when the fraction is expanded as a series in a it is identical to the raw series. This approach markedly improves the 
results for p. For example, for the random walk, the Pade of the 14th order series appears to better the 25th order 
raw series obtained from the matrix method (both raw series agree, of course, up to 14th order). However, in order 
to get accurate continuum results for 9 we add 1 further term to the Pade (either numerator or denominator) whose 
coefficient is chosen so that the exact constraint p{\) — 1 is satisfied. This serves to give reasonably accurate estimates 
of the continuum persistence exponent. For example, for the random acceleration problem we find 9 = 0.2506(5) from 
the Pade approach, compared to the exact result of 1/4. 

For certain sufficiently smooth processes the derivative of 9d{XT) at XT — is zero [24]. We can thus add a further 
term to the Pade to impose this constraint, and markedly improve the accuracy of our estimate of the continuum 9. 
The diffusion equation in all dimensions and the d n x/dt n = r](t) processes for n > 3 are all suitable. Table VI shows 
the continuum results for diffusion in 1-3 dimensions as reported in [25], with numerical results and also the singly 
constrained and IIA results for comparison. For d — 1 the IIA is slightly better, but the correlator expansion is more 
accurate for d = 2,3. Furthermore, we obtain estimates of the errors and, by going to higher order we may improve 
our results. Table VI shows the continuum results for the d n x/dt n = r/(t) processes with 3 < n < 10. Figure 9 shows 
how 9 n varies with n. Notice in particular that 6 n — 6^ oc 1/n for n > 20, and that 9^ is identical to that of 2-d 
diffusion (since the correlators are identical). 





PadclCR 


Padc2CR 


numerical 


IIA 


X, 

1- d diff 

2- d diff 

3- d diff 


0.2506(5) 
0.119(1) 
0.187(1) 
0.24(3) 


0.1201(3) 
0.1875(1) 
0.237(1) 


1/4 (exact) 
0.12050(5) 
0.1875(1) 
0.2382(1) 


0.2647 
0.1203 
0.1862 
0.2358 



Table VI. Results for the continuum persistence exponent 9 for the random acceleration problem and the diffusion 
equation in 1-3 dimensions. PadelCR is the Paded results with 1 constraint, Pade2CR has 2 constraints. Numerical 
[13] and IIA results are shown for comparison. 
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0.19864(3) 


0.19813 


8 


0.19707(3) 


0.19642 


9 


0.19589(3) 


0.19514 


10 


0.19496(3) 


0.19414 



Table VI. Results for n against n for small n. The Paded correlator expansion with 2 constraints is shown along 
with the Independent Interval Approximation. The Paded results are an average of suitable Pades of order 14 to 10. 
Note that for n = 2 (the random acceleration problem), the IIA gives 9 = 0.2647 whilst the analytical result is 1/4. 
For n — > oo (the diffusion equation), the Paded result with 2 constraints is 0.1875(1), the numerical result is 0.1875(1) 
and the IIA result is 0.1862. 
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1/n 
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n 

FIG. 9. Plot of n against n for small n. The n = 1, 2 
results are omitted as they are known analytically to be 1/2 
and 1/4 respectively and only one constraint may be imposed 
on the Pade in these cases. Note that 9 n goes to the contin- 
uum result of 0.1875(1) (solid line) rather slowly, in fact as 
1/n, see inset. The results were obtained using Pades with 2 
constraints, an average being taken of suitable Pades of or- 
der 14 to 10. Inset: Plot of 6 n against 1/n showing that 
n — Qoo oc 1/n for n > 20. 

And so, by the use of these two constraints we have been able to extend a series expansion about AT = oo all the 
way to the AT — > limit. The ability to do this does however depend on the correlator. First, if the process is 'rough', 
i.e. 1 — C(T) oc T 13 + . . . with < (i < 2, so that the probability distribution of the time T between two zero-crossings 
behaves as P X {T) oc T a + . . ., with a < 0, then we have shown [24] that 0(0) - 6»(AT) - Z^T 1 + Q for small T. For 
the random walk, for example, a = —1/2 and so we get a square root cusp in #£>(AT) for AT — > which the series 
expansion about AT = oo cannot reproduce. Consider processes such as 

r)h 

-^ = -(-v 2 y /2 h+ v (t) (60) 

describing linear interface growth, where ry(x, t) is delta correlated in space and time. The normalized autocorrelation 
function of h(x,t) has, for d < z, the form 

C(T) = cosh(T/2) 2/3 - | sinh(T/2)| 2/3 = 1 - |T/2| 2/3 + . . . (61) 

where T = lnt as usual, (i — (1/2) (1 — d/z) and d is the spatial dimension [28]. Since (3 < 1/2 this process will always 
be 'rough', so extrapolation of the scries to the continuum limit is not possible. Secondly, the correlator may not be 
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easily expandable in some suitable variable such as the e _Ar / 2 used earlier. Another process, fractional Brownian 
motion, defined as a Gaussian process x(t) with stationary increments ([x(t) — x{t')] 2 ) oc \t — t'| 2/3 , has normalized 
correlator [29] 

C(T) = cosh(/3T) - i 

where T = lnt as usual. For general (3 there are two incommensurate variables e~^^ and e _Ar / 2 , making it difficult 
to construct a controlled expansion. 

Finally, when applying the 2 constraints to the 3-dimensional diffusion problem, we have been unable to numerically 
solve the 44 simultaneous nonlinear equations required to construct the expansion to 0(a 42 ). Thus we have only gone 
up to a 29 in this problem. This is not however an insuperable difficulty. 

Note however that, even when we cannot get continuum results, for AT large the expansion will always work, as 
just substituting in the raw correlator is good enough. 

Having introduced the correlator expansion for the calculation of persistence exponents, in the following sections 
we will extend it to calculate properties of the occupation-time and crossing-number distributions. 

VII. OCCUPATION-TIME DISTRIBUTION 

The occupation-time distribution, considered for the continuous case in [9-11], is the probability distribution 
R(t(T)), where 

T(T) = ±£dT'G(X(T')) (63) 

and <d(X) is the Heaviside step function. For a symmetric distribution of zero mean, Rt(t) is symmetric about 
t = 1/2. Then i?r(0) and -Rt(I) give the persistence probability P(T) introduced earlier. The discrete-sampling 
equivalent, R n , s , is the probability that X(T) has been found to be positive at exactly s out of the n samplings. Thus, 

n 

s(n)/n = r(n) = -^OLY^AT)] . (64) 
n i=i 

Writing Q[X(iXT)] = (1 + cr 4 )/2, where a { = signLY(iAT)], we get 

(r(n)) = \ (65) 

and 

1 n n 

< r(n)2 > = 4 + W £ £ sin_1 [ C ^ - Ml ( 66 ) 

where we have used the result that {<Ji<jj) — (2/7r)sin _1 [C(|z — j|AT)]. If we choose as before to work to a given order 
in the correlator, we need only evaluate the sum up to that order. Taking the large n limit, we can change the sum 
to: 

(r{nf) = \ + ^ + ^±^[C{kKT)] (67) 
fc=i 

where o is the order to which we wish to work. It has been pointed out [16] that for n large, the widely separated 
(in time) parts of the time series become uncorrelated and, following the central limit theorem, R n _ s is Gaussian for 
s close to 1/2 with standard deviation given by equation (67). We will use this as a check of our final result for R n , s - 
To find R n>s we must sum over all 'paths' involving s positive samplings (and n—s negative ones). So the probability, 
Rn tS to find s positive values from n samplings is 

Rn, S = (S 2s _ niE ai ) = (^ £ ft 1 + ( 6g ) 

\ {£i=±l} * »=1 / 



2 sinli | — 



(62) 
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where 5 a ^p is the Kronecker delta function which we choose to write in analytic form as a Cauchy integral 

1 f dz 



<W = ^ f ^=^, (69) 



where the integration contour encircles the origin. Substituting this into eqn. (68) gives, 



\ {ei=±l}*=l / 



Summing over the e^s gives: 

1 f dz 



R n 



2-rti T z 2s - n+1 




The term which is averaged over is identical to that of the normal persistence calculation apart from the factor 
(z 2 — l)/(z 2 + 1) associated with each Oi. Making this minor change to the previous calculation of persistence gives 
us a term T n (z) where before we had p n , and so T(oo) = p. Replacing s by rn where < r < 1 and anticipating that 
Rn,s ~ [p( r )] n f° r n large gives us an expression for p(r) which we can evaluate by steepest descents: 



WW = ^ / f ex P [n(ln(l + 1) - rlni + InT(t))] 



(72) 



where we have replaced z 2 by t, and T(t) = Y(\/i). As a simple check, at zeroth order T is 1/2 and the method of 
steepest descents gives a saddle-point value 4°^ = r/(l — r), and 

Rn,s ~ ^exp{-n[rlnr + (l-r)ln(l-r)]}. (73) 
This is the same as the combinatorial result, 

r --=2K;i) <74) 

when expanded to leading order for large n using Stirling's formula. Note that there is hence also a \fn term in R n , s 
which we ignore relative to the exponential for n — > oo. 

We use the method of steepest descents in the following way. Having found the position of the saddle point to 
zeroth order, we substitute it into the right-hand side of the general saddle-point equation 

t s = j-~ ts ^ + ts) T'(t s ), (75) 

1 — r 1 — r 

where T'(i) = dT/dt, and thus find t s to first order, and so on recursively up to 10th order. Substituting t s into the 
exponent of equation (72) gives an analytic expression for R ntS in the large n limit and hence p(r). Just as in the 
expression for persistence, the expression for p(r) is rather long and it was only possible to find p(r) analytically to 
10th order. 

As stated in the previous section, for r close to < r >= 1/2, p(r) approximates to a gaussian distribution, 



p(r) cx exp 



1 [r-{r)f 
2n (r 2 ) - <r\ 



(76) 



Thus one expects that the quantity lim r ^( r ) (r— (r)) 2 /[— 2n\np(r)] should equal the variance of r calculated previously, 
and indeed these two quantities agree term by term to 10th order, providing a useful cross-check. 

We apply our result to the random walk, random acceleration, and diffusion from random initial conditions in 1-3 
dimensions. We recall that R n s is the probability for n measurements of X to return s positive values. We have 
shown that for n — > oo, s — > oo with r = s/n fixed it has the form R n ^ s ~ [p(r)] n , which can be written in the 
alternative form R nyS ~ exp[—0£>(r)T], where T = nAT as usual and 0d = — lnp(r)/AT. 
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Plots of 6o(r) against r for various values of AT are shown in figures 11, 12, 13, 14 and 15. For the diffusion 
equation we are able to Pade the series and apply 2 constraints, giving us good estimates for the continuum 9(r), 
i.e. the limiting value of 9o(r) as AT — > 0. Plots of 9(r) are also shown in figure 10. The second constraint, that 
d#(r)/c?AT|AT=o = for sufficiently smooth processes (including diffusion), comes from a similar argument to that 
given earlier [24] for standard persistence: as AT is increased from zero, the first correction to 9 comes from the 
contribution of a path that is the same as a contributing path in the continuum, apart from one undetected double 
crossing which (to lowest order in AT) gives a correction to 9 of order AT 2 and thus d9(r) / 'cZAT|at=o = 0. 

The function 9{r) is the large-deviation function for the occupation-time distribution. Close to r = (r) = 1/2, it is 
quadratic in r — (r) . The probability distribution P r (r) of r is given by P r (r) oc [p{r)] n = cxp[— (l/2)(r — (r)) 2 /((r 2 ) — 

(r) 2 ) for r near (r). This means that the typical fluctuations in r around the mean are of order n -1 / 2 for large n since 
the variance is proportional to 1/n. The full function 9{r) is required to determine the probability of large deviations 
from the mean, where the fluctuations are non-Gaussian. 

We end this section by noting that the full large deviation function #o(r) associated with the occupation-time 
distribution was computed analytically [18] for the intrinsically discrete process 

ipi = cos(cu)fa + sm(u)fa-i , (77) 

where the fa are independently distributed gaussian random variables. This process appears as a limiting case of 
the diffusion equation on a hierarchical lattice [30] and also appears in the one dimensional Ising spin glass problem 
[18,31]. Exact results were obtained for the case lo = tt/A. Interestingly, these results turn out to be independent 
of the distribution of fa provided that it is symmetric. We have also obtained the the large deviation function for 
w = 7r/4 by the correlator method. The comparison with the exact results is shown in figure 16. 

0.25 | . 1 . 1 . 1 . 1 . 1 
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FIG. 10. Plot of the continuum large deviation function 
6(r) against r for the diffusion equation in 1, 2 and 3 dimen- 
sions (bottom to top respectively). 6(r) is symmetric about 
r — 1/2. The results were obtained using Pades with 2 con- 
straints, an average being taken of suitable Pades of order 10 
to 7 for 1 and 2 dimensions and of order 7 to 6 for 3 dimen- 
sions. 
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FIG. 11. Plot of 9d(t) against r for the random walk with 
cxp(-AT/2) = 1/2, 1/4, 1/8, 1/16, 1/256 and 1/2 40 (top to 
bottom respectively). The curves were plotted from the raw 
series in exp(— AT/2) to 10th order. Note that r = 0, 1 cor- 
responds to ordinary discrete persistence and that the curves 
are symmetric about r = 1/2. 
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FIG. 12. Plot of 9d{t) against r for the random acceler- 
ation with exp(-Ar/2) = 1/2, 1/4, 1/8, 1/16, 1/256 and 
1/2 40 (top to bottom respectively). The curves were plotted 
from the raw series in exp(— AT/2) to 10th order. Note that 
r = 0, 1 corresponds to ordinary discrete persistence and that 
the curves are symmetric about r = 1/2. 
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FIG. 13. Plot of 9o{r) against r for diffusion in 1 dimension 
with exp(-AT/2) = 1/4, 1/8, 1/16, 1/256 and 1/2 40 (top to 
bottom respectively). The curves were plotted from the raw 
series in exp(-AT/2) to 10th order. The exp(-AT/2) = 1/2 
curve is not shown as the raw series had not converged at 10th 
order. Also shown are the results for the continuum limit from 
figure 10. Note that r = 0, 1 corresponds to ordinary discrete 
persistence and that the curves are symmetric about r = 1/2. 
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FIG. 14. Plot of 0d{t) against r for diffusion in 2 dimension 
with exp(-Ar/2) = 1/2, 1/4, 1/8, 1/16, 1/256 and 1/2 40 
(top to bottom respectively). The curves were plotted from 
the raw series in exp(— AT/2) to 10th order. Also shown are 
the results for the continuum limit from figure 10. Note that 
r = 0, 1 corresponds to ordinary discrete persistence and that 
the curves are symmetric about r = 1/2. 
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FIG. 15. Plot ol 0d(t) against r for diffusion in 3 dimension 
with exp(-AT/2) = 1/4, 1/8, 1/16, 1/256 and 1/2 40 (top to 
bottom respectively). The curves were plotted from the raw 
series in exp(-3AT/2) to 10th order. The exp(-AT/2) = 1/2 
curve is not shown as the raw series had not converged at 10th 
order. Also shown are the results for the continuum limit from 
figure 10. Note that r = 0,1 corresponds to ordinary discrete 
persistence and that the curves are symmetric about r = 1/2. 




FIG. 16. Plot of 6o(r) against r for the intrinsically dis- 
crete process ipi = (<f>i + <^_i)/\/2 where the 4>iS are indepen- 
dent identically distributed symmetric random variables. The 
curves show the raw correlator result to 10th order (dotted 
line) and the exact result [18] (solid line). The curves differ by 
a maximum of 0.00575. Note that r = 0, 1 corresponds to or- 
dinary discrete persistence and that the curves are symmetric 
about r = 1/2. 
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VIII. OCCUPATION-TIME PARTIAL SURVIVAL 



Here we consider the discrete occupation-time partial-survival probability, R n {p). Let us suppose that the process 
'dies' with probability 1—p whenever X is sampled to be positive. Then R n (p) is defined to be the probability of the 
process surviving n samplings if the variable Xi survives being sampled as positive with probability p. Samplings as 
negative are always survived. Thus, 

n 

Rn(P) = Y,P SR ^ ( ?8 ) 

and so R n (p) is also the generating function for i?„ jS since 



= 7\^ p=Mp) (79) 



or alternatively 

Rn,s = j ^Rn{p) (80) 

Also, writing p s as exp(slnp) and expanding the exponentials gives: 

w(p) = E ¥ V) e (81) 



where < s J > c is the jth cumulant of the occupation time, s, and we have used R n (p) — [p(p)] n , which, as for 
persistence (but unlike R n , s ) is true for any n provided that n is larger than the largest diagram involved in the 
evaluation of R n {p)- Thus calculating R n {p) gives us another method for finding the moments of the number of 
crossings and also although the evaluation of by the contour integration is entirely equivalent to the 

previous method and differentiating R n {p) s (= rn) times becomes unfeasible for large n. 
R n {p) is found in a similar way to before, by summing over all possible 'paths', 

\ e»=±li=l v / 

where Uj = sgn(X i ) as usual, and the average is over the variables Xi (i = 1, . . . , n). Thus we get 

^(p)=(^(p+ i ) B n( i +fri^)) ( g3 ) 

which is the same as the calculation for normal persistence except that we include a factor (p— l)/(p+ 1) with each at 
and and overall factor (p+ 1)™. Thus we can find p(p) — cxp[— 9(p)] to order 14. This is done and results for diffusion 
in 1-3 dimensions are shown in figures (17, 18, and 19). Note that, for /XT — > 0, a positive excursion by the underlying 
continuous process will survive with zero probability since the number of samplings — > oo. Thus 0(p)|at=o is just the 
continuum persistence exponent. It is therefore possible to improve 9{p) for XT small by applying 0(p)|at=o = as 
a constraint, in addition to the standard constraint p(p)\xr=o = 1- 
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FIG. 17. Plot of 6 D (p) against exp(-AT/2) for diffusion in 
1 dimension with p = 0, 1/4, 1/2, 5/8, 3/4 and 7/8 (from 
the top respectively). Qd{p) has been constrained to give the 
persistence result at the continuum. The curves are produced 
from averages of suitable constrained Pades, although in prac- 
tice the various Pades are indistinguishable. 
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FIG. 18. Plot of 6 D (p) against exp(-AT/2) for diffusion in 
2 dimensions with p = 0, 1/4, 1/2, 5/8, 3/4 and 7/8 (from 
the top respectively). Qd{p) has been constrained to give the 
persistence result at the continuum. The curves are produced 
from averages of suitable constrained Pades, although in prac- 
tice the various Pades are indistinguishable. 
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FIG. 19. Plot of 9 D (p) against exp(-AT/2) for diffusion in 
3 dimensions with p = 0, 1/4, 1/2, 5/8, 3/4 and 7/8 (from 
the top respectively). 9d(p) has been constrained to give the 
persistence result at the continuum. The curves are produced 
from averages of suitable constrained Pades, although in prac- 
tice the various Pades are indistinguishable. 



A further check is provided by using p(p) to generate the first two cumulants. The results agree term by term to 
10th order with the method used to calculate the mean (trivially) and the variance. 

We next compare the results obtained by the correlator method to an exactly solvable case, namely the discrete 
process in Eq. (77) for u) = w/4. In this case, an exact expression of the exponent 8(p) is known [18], 



9(p) = 



(1-P) 



(84) 



A comparison of this exact result with the one obtained by the correlator method is shown in figure 20. 




0.6 
P 

FIG. 20. Plot of 9d{p) against p for the process 
ipi = {4>i + 0i-i)/v / 2- The exact (solid line) and raw correla- 
tor expansion (dotted line) results are shown, the two curves 
being indistinguishable except near p = 
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In the last 2 sections we have examined the occupation-time statistics. The occupation time depends only on the 
signs of X(i/XT) at each i, that is, it is local. This meant that we merely had to attach additional factors to each 
local X{i/XT). In the next sections we will be studying the number of crossings, so we must consider the signs of 
both X(iAT) and X((i + 1)AT). Thus the problem is not local in the sense used above and we cannot just attach 
additional factors to each X(iAT). The solution, as explained in the next section, is to attach additional factors to 
the lines connecting two Xs in the diagrammatic notation. 



IX. DISTRIBUTION OF CROSSINGS 



We now apply the correlator expansion to calculate the distribution of crossings of an arbitrary GSP. We start from 
the calculation of the persistence. The method is the same up until we assign factors to the diagrams on the lattice. 
We wish to calculate the probability of m detected crossings in n samplings, P„ >m , rather than just the probability 
of no crossings which was the persistence calculation. To do this, we sum over all the possible ways in which those 
to crossings could occur. Furthermore, we note that if we have a line in a diagram connecting two vertices, and s 
crossings occur between these two vertices, then the factor C(jAT) associated with it from the persistence calculation 
should also have a factor (— l) s associated with it. Consider the diagram shown in figure 21. 



FIG. 21. A 4th order diagram. 

Besides the enumerations done for the persistence calculation we must consider the following four cases. 
1: No crossings occur on the sites where the diagram is placed. This would occur with probability 

(n-ro - l)(rc-ro-2) 

(n-l)(ra-2) (85j 

and there are no sign changes. 

2: There is one crossing between the first and second vertices. This occurs with probability 

777(77 — m — 1) 

(77 - 1X77 - 2) (86) 

and there is a factor (— l) 3 associated with it. 

3: There is one crossing between the third and fourth vertices. The probability of this occurring is as above and there 
is a factor (—1) associated with it. 

4: There is one crossing between the first and second vertices and one crossing between the third and fourth vertices. 
This occurs with probability 



777(777 — 1) 
(77 - 1)(t7 - 2) 

and there is a factor (— 1) 3 (— 1) associated with it. Hence this diagram has an additional factor 

(77 — 777 — l)(n — 777 — 2) — 2m(n — 777 — 1) + 777(777 — 1) 
(t7-1)(t7-2) 

over and above that for the persistence calculation. Furthermore, there is an overall factor 

77-I 
m 



(87) 



(88) 



(89) 



accounting for all the ways in which the m crossings can occur on the whole lattice. Note that, although we use 
the term 'probability', when run over the whole lattice (multiplied by the factors (77 — 2) (77 — 3)/2!) each probability 
becomes the exact number of ways in which the corresponding event occurs. Thus by introducing these extra rules 
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when enumerating the diagrams on the lattice we are able to calculate P n 
crossings occurring in n samplings. For n large we expect that, as usual, 



p 

1 n,'. 



the probability of exactly m detected 



(90) 



and so we find p m as 



,. n+l,m+m/n 

= hm — — . 



(91) 



This has been done, although due to the additional factors it was possible only to go to 10th order due to memory 
constraints. It has been checked that the result agrees term by term with the normal persistence calculation for m = 
and with alternating persistence for m — n. Figures 22, 23, 24, 25, and 26 show p(r) against r where rn = m for 
various values of exp(— AT/2) for the random walk, random acceleration, and diffusion in 1-3 dimensions. Notice that, 
for r = (r) = 1/2 — (1/n) sin -1 C(AT), p(r) = 1. Close to this point, p(r) approximates to a Gaussian distribution, 



p(r) oc e 2(<^>-<,-> 



(r—<T»* 

~T: 



(92) 



where the variance < r 2 > — < r > 2 agrees term-by-term with the calculation in the following section. Remember 
that there is a next-to-leading term (pre-exponential factor), sqrtn in -P„, m , as can be seen from considering the 
AT = oo (lowest order) case: 



P A. 



n- 1 

m 



1 



1 



2 (1 - r ) 1 " r r 7 



(93) 



Note that we are considering the number of detected crossings per sampling (0 or 1). As AT — > the fraction of 
(detected) crossings will go to zero. Because of this, 6{r) — > oo for AT — > except for the r = case which is just 
standard persistence. As always, < p(r) < 1, and we choose to plot p(r) rather than 9(r). Wc have not applied any 
constraints to the series. 
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FIG. 22. Plot of p(r) against r for the random walk with 
exp(-AT/2) = to 8/10 in steps of 1/10 (top right to bottom 
right respectively). Note that p(r) is 1 at the mean value of 
r given by (r) = 1/2 - sin" 1 [G(AT)]/7r. The plots are of the 
raw series. 
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FIG. 23. Plot of p(r) against r for the random acceleration 
process with cxp(-Ar/2) = to 7/10 in steps of 1/10 (top 
right to bottom right respectively). Note that p(r) is 1 at the 
mean value of r given by (r) = 1/2 — sin _1 [C(Ar)]/7T. The 
plots are of the raw series. 
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FIG. 24. Plot of p(r) against r for diffusion in 1 dimension 
with exp(-Ar/2) = to 6/10 in steps of 1/10 (top right to 
bottom right respectively). Note that p(r) is 1 at the mean 
value of r given by (r) = 1/2 - sin -1 [C(AT)]/ir. The plots 
are of the raw series. 
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FIG. 25. Plot of p(r) against r for diffusion in 2 dimensions 
with exp(-AT/2) = to 7/10 in steps of 1/10 (top right to 
bottom right respectively). Note that p(r) is 1 at the mean 
value of r given by (r) = 1/2 - sin -1 [C(AT)]/7r. The plots 
are of the raw series. 
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FIG. 26. Plot of p(r) against r for diffusion in 3 dimensions 
with cxp(-AT/2) = to 3/10 in steps of 1/10 (top right to 
bottom right respectively). Note that p(r) is 1 at the mean 
value of r given by (r) = 1/2 - sin" 1 [C(AT)]/7T. The plots 
are of the raw series. 

Recently one of us [5] calculated the distribution of crossings and partial-survival probability of the intrinsically 
discrete process tpi given by Eq. (77) for the special case w — 7r/4. The correlator of the process is 

C{i — j) = 5i,j + cosu) smw(8ij-i + Sij+i). (94) 

Substituting this into the correlator expansion gives the result shown in figure 27 for comparison with the analytic 
result. The agreement is good for r small but for r > 0.73 the series has not yet converged by 10th order. This shows 
up in the way that p(r) changes as the order is increased from 1 to 10. For r small there is oscillatory convergence 
whilst for r large the convergence is monotonic or, for r > 0.73, has not occurred. The fact that convergence does 
not occur for r large is presumably because the series is less good for large numbers of crossings. This also occurs, for 
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example, for the random acceleration problem where the alternating persistence (r = 1) result converges more slowly 
than the standard persistence (r = 0) result. That it fails so badly whilst the small r result is acceptable is surprising. 
Nevertheless, by checking whether or not the series converges as the order is increased to 10, we can tell whether the 
result is reliable. For the case ui = n/12, for which C(i — j) = 5ij + j(Sij-i + <5 ij+ i), the series has converged for all 
r although there is no analytic result for this case. In fact the case studied is the one for which the correlator takes its 
largest possible value. Also notice that the p(l) = result is due to the requirement that, in order that tpi alternate 
in sign, the magnitude of <f>i must increase every time step. Thus P n , n decays as 2~ n /n\, which is faster than a power 
of n, implying p(l) = 0. For other values of the coefficients of and <$H-\ this is not the case and presumably p(l) 
is non-zero. 




r 



FIG. 27. Plot of p(r) against r for the process 
ipi = {4>i + 0i-i)/v / 2 where the (f>iS are independent identi- 
cally distributed symmetric random variables. The solid line 
is the numerical evaluation of the exact result and the dashed 
line is the result of the correlator expansion. The agreement 
is good until r ~ 0.73, and becomes badly wrong as r — > 1 
(see text). 

In this section we have calculated the distribution of crossings, P„, m , to 10th order in the correlator by extending 
the diagrammatic technique. In the next two sections we derive the standard deviation of the number of crossings 
and then use P n , m to calculate F n (p), the partial survival of crossings probability. 

X. THE VARIANCE OF THE NUMBER OF CROSSINGS 

The number of detected crossings in n samplings, m, is (up to an end effect that is negligible for large n), 

n 1 

m = rn = ^2-(l-aia i+1 ) (95) 

where Oi = sign{X (i&T)]. This gives 

(r) = i-isin- 1 [C(AT)] (96) 
as derived in [26]. One may further attempt to evaluate the variance of m, 

_^ n n 

cr 2 /n = ((m 2 ) - (to) 2 ) jn = — ( <J i C7 i+i <J J <7 J+i> - (<W+i) (97) 

j=i j=i 
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which involves calculating the connected 4-vertex diagrams in the correlator expansion. The calculation is essentially 
identical to that of section VI apart from the extra cases of i = j and i = j ± 1 and the result to 14th order may be 
read off. Figures (22 ,23 ,24 ,25 ,26 ,27) show p(r) against r for various processes and values of AT. Close to r = (r), 
p{r) is given by 



p(r) ~ exp 



1 (r-(r)) 2 
2n ( r 2) - (rf 



(98) 



and comparison of (r— (r) ) 2 /[— 2nlnp(r)] agrees term by term to 14th order in the correlator with the direct calculation 
of (r 2 ) — (r) 2 , providing a useful cross-check. 

The result for (r 2 ) — (r) 2 for the random walk also agrees with that of equation (20) (the matrix method). The 
variance for various processes is shown in figure 28. 
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FIG. 28. Plot of a 2 against exp(-AT/2). The curves are, 
from the top right; the linear growth equation (60), the ran- 
dom walk, random acceleration, and diffusion in 3, 2 and 
1 dimensions. The curves are the raw series in powers of 
exp(— AT/2) to 14th order in the correlator and are plotted 
only as far as their series have converged. 

Thus we have found the variance of the number of crossings for an arbitrary process to 14th order in the correlator. 
This involved calculating the 4-vertex diagrams only, and therefore it is entirely feasible to go to higher orders since 
the 4-vertex diagrams are relatively simple. Notice also that the variance only contains even orders, as one would 
expect from the 4-vertex diagrams. 

In the following section we complete our calculations by finding the partial-survival probability for an arbitrary 
GSP, this also being the moment generating function. The results will be shown to agree with those of the current 
section. 

XI. DISTRIBUTION OF CROSSINGS PARTIAL SURVIVAL 

As for the specific case of the random walker (Section V), we may consider the partial-survival probability F n (p), 
the probability of surviving up to the nth sampling if each detected crossing is survived with probability p. As stated 
in Section V, this is also the generating function for P n , m and the cumulant generating function [4]: 



F n ( P )= J2p mP n, 



(99) 



m=0 



and 
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^(P) = ET<» J >c ( 10 °) 
3=1 ■'' 

where (r 7 ') is the jth cumulant. From equation (99) it can be seen that F n (p) is a sum of terms containing (^)m s p m 
where s is some positive integer. These can be simply evaluated to give an expression for F n (p) and hence p p . As for 
the occupation partial survival, F n (p) = p{p) n = exp[— 9(p)n] for all n even though this is not true for P n , m - 

For rough processes, the continuum partial survival is the same as persistence, since any crossing entails an infinite 
number of crossings and thus non-survival. For smooth processes however, calculation of p(p) and hence that of 
8{p) = —ln[p(p)} is non-trivial. For the random acceleration problem the exact result is [32], 

^)=i( i -^ sin "i)- ( io1 ) 

Whilst for the intrinsically discrete process (equation (77)) it is [5], 

^'"P"^ 2 )' (102) 

Figure 29 shows this result and the raw series result for comparison. 

For general processes we apply the constraint to the series that p(p)\&r=o = 1- For sufficiently smooth processes, 
as before, the first correction to 9(p) near AT = will be of order 1ST 2 . We apply this additional constraint to the 
appropriate processes and the corresponding continuum results are shown in figures 30, 31. 
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FIG. 29. Plot of 8(p) against p for the process 
ipi = (<f>i + 0i-i)/v / 2- The solid line is the exact result and 
the dashed line is the result of the correlator expansion. They 
differ by a maximum of 0.00576 (at p = 0). 
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FIG. 30. Plot of 6{p) against p for (top to bottom) the 
d n x/dt n = r)(t) process with n — 2 (random acceleration), 
n = 3, n = 4, and n — > oo (equivalent to diffusion in 2 
dimensions). For n — 2 (top curves) the exact result (101) is 
shown (dashed line) along with the IIA (dotted) and the Pade 
with 1 constraint (solid line). For the other cases, the IIA 
(dotted) and Pade with 2 constraints (solid line) are shown. 
The correlator results are an average of suitable Pades of order 
10 to 7. Note that the IIA is rather inaccurate for the random 
acceleration, but improves as n increases. 
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FIG. 31. Plot of 8(p) against p for diffusion in 1, 2 and 3 
dimensions (bottom to top). The IIA results (dotted) and 
the correlator results (solid lines) are plotted. The correlator 
results are an average of suitable Pades of order 10 to 7 for 1 
and 2 dimensions and 6 to 5 for 3 dimensions. 

The results from the matrix method partial survival for the random walk and random acceleration agree with those 
of the correlator expansion term by term to within the numerical error of the matrix method. Also, using F n (p) as a 
generating function, we find that the first two cumulants agree term by term to 10th order with the results of section 
X. We are also able to calculate higher cumulants. We have also used the method of steepest descents to calculate 
p(r) from p(p) as a further cross-check. 
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This evaluation of the crossing partial survival completes our calculations. 



XII. CONCLUSION 

In this paper we have used the correlator expansion to calculate several properties of an arbitrary discrete or 
discretely sampled Gaussian Stationary Process. The expansion in powers of the correlator works well when the 
variables are relatively weakly correlated. For stronger correlations the series expansion does not converge. We are 
however able, for the case of an underlying continuous and sufficiently smooth process, to extrapolate our results all 
the way to the continuum by using the Pade Approximant with two constraints. Thus even in the continuum we have 
been able to calculate the persistence exponents, the occupation-time exponents and the partial survival of crossings 
exponents to high precision. These results compare well with those of the Independent Interval Approximation, the 
other general method. In most cases they are more accurate, however, and they also give an estimate of the error of 
the result which the IIA does not. Furthermore, by calculating higher orders the results may be improved. We believe 
we have demonstrated convincingly that the correlator expansion is the method of choice for calculating persistence 
properties of Gaussian stationary processes. 
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